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ABSTRACT 

This  technical  memorandum  describes  the  fortran  computer  program  WADAl  (WAve 
Diffraction  Analysis — version  1),  which  calculates  the  pressure  induced  by  the  diffraction  of 
regular  incident  waves  on  arbitrary-shaped  two-dimensional  cylinders  fixed  in  the  free  surface 
of  deep  water.  Based  on  the  source-distribution  method,  WADAl  solves  the  integral  equation 
for  the  linear  diffraction  problem  by  Frank's  close-fit  method. 


RfiSUMfi 

Le  present  rapport  technique  decrit  le  programme  fortran  WADAl  {analyse  de  la  diffrac¬ 
tion  des  vagues,  premiere  version)  qui  permet  de  calculer  la  pression  induite  par  la  diffraction 
de  vagues  incidentes  regulieres  sur  des  cylindres  d  deux  dimensions  et  de  formes  arbitraires 
fixes  a  la  surface  lihre  d’eaux  profondes.  En  faisant  appel  d  la  methode  de  distribution  des 
sources,  le  programme  WADAl  resout  I ’equation  integrate  du  probleme  de  diffraction  lineaire 
au  moyen  de  la  methode  de  I’ajustement  fin  de  Frank. 
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amplitude  of  incident  wave 
breadth  of  cylinder  at  waterline 
Green  function 
gravitational  acceleration 
influence  coefficients 
imaginary  part  of  a  complex  variable 

integral  of  dG/dn  along  Sj  for  field  point  located  at  the  node 

dummy  variable  for  integration 

normal  unit  vector  to  S  pointing  out  of  fluid 

normal  unit  vector  to  5,;  n,  = 

X-  and  y-components  of  n, 

time-dependent  incident-wave  pressure  (=p;e~*“*) 
time-dependent  diffracted-wave  pressure  (= 
time-dependent  total  wave-induced  pressure  (=  Pi  +  Pd) 

Cauchy  principal  value 
complex  amplitude  of  incident-wave  pressure 
complex  amplitude  of  diffracted-wave  pressure 
source  strength 

source  strength  for  the  segment  Sj 
real  part  of  a  complex  variable 
cylinder  contour  below  calm-water  line 
segment  of  5 
length  variable  along  S 
draft  of  cylinder 
time 

normal  velocity  of  incident  wave  at  Si 
coordinates  of  a  field  point 
complex  variable  (=  i  -1-  »y),  or  a  coordinate 
wave  direction  relative  to  the  z-axis 
Kronecker  delta  function 

integral  of  G  along  Sj  for  field  point  located  at  the  node 
complex  variable  (=  C  +  ir}) 


* 


f  complex  conjugate  of  f 

Tj  coordinates  of  a  source  point 

frequency  parameter  (=  i>Bf2) 

V  wave  number  (= 

p  mass  density  of  fluid 

time-dependent  incident-wave  potential  (= 

$£)  time-dependent  difFracted-wave  potential  (= 

time-dependent  total  potential  (=  $/-!-$£)) 

4>  phase  angle  between  wave-induced  pressure  and  incident  wave 

(pj  complex  zjnplitude  of  incident-wave  potential 

<Pb  complex  amplitude  of  difFracted-wave  potential 

u  wave  frequency 

V  gradient  operator 


Superscript 

(m)  mode  of  pressure  excitation 

m  =  1  horizontal  excitation  (eisymmetric  about  the  y  axis) 
m  =  2  vertical  excitation  (symmetric  about  the  y  axis) 

Subscripts 

i,j  indices  for  contour  segments 

n  indices  for  irregular  frequencies 


I 


V 


1.  INTRODUCTION 


This  technical  memorandum  describes  the  fortran  computer  program  WADAl  (WAve 
Diffraction  Analysis — version  1)  developed  as  part  of  DREA’s  on-going  research  to  analyze 
the  interactions  of  ocean  waves  with  marine  vehicles.  WADAl  calculates  the  hydrodynamic 
pressure  induced  by  an  incident  sinusoidal  wave  and  diffracted  wave  on  arbitrary-shaped  two- 
dimensional  cylinders  fixed  in  the  free  surface  of  deep  water. 

WADAl  solves  the  linear  diffraction  problem  by  the  source-distribution  method.  It  uses 
the  Green  function  that  is  a  particular  fundamental  solution  of  the  two-dimensional  Laplace 
equation  and  directly  satisfies  the  boundary  conditions  on  the  free  surface,  on  the  sea  bottom, 
and  at  infinity.  The  amount  of  required  input  data  is  minimal  since  only  the  boundary  condition 
on  the  body  surface  remains  to  be  satisfied;  however,  the  computation  of  the  Green  function 
and  its  gradient  is  complicated.  Also  the  method  fails  to  give  a  unique  solution  at  “irregular 
frequencies,”  but  these  usually  occur  outside  the  range  of  frequency  of  practical  interest. 

The  mathematical  formulation  of  the  linear  diffraction  problem  is  described  in  Section  2. 
The  solution  technique  adopted  in  WADAl,  a  simple  collocation  method  originally  used  by 
Frank  is  outlined  in  Section  3.  Comments  on  the  “irregular  frequencies”  are  made  in  Section 
4.  Descriptions  of  the  input  and  output  are  given  in  Sections  5  and  6,  respectively.  Concluding 
remarks  are  given  Section  7.  Sample  input  and  output  are  included  in  Appendices. 

2.  STATEMENT  OF  THE  MATHEMATICAL  PROBLEM 

The  water  is  assumed  to  be  deep  and  inviscid,  and  its  motion  irrotational.  The  wave 
amplitude  is  assumed  to  be  small  compared  with  a  characteristic  length  of  the  body,  say,  the 
half-breadth  of  the  cylinder,  so  that  linearization  is  valid.  In  Fig.  1,  the  Cartesian  coordinates 
(i,y,  z)  are  fixed  to  the  cylinder,  which  is  assumed  to  be  of  infinite  length  in  the  z-direction 
with  a  uniform  cross  section.  The  iz-plane  lies  in  the  calm-water  surface.  The  y-axis  points 
vertically  upwards.  The  origin  is  located  at  the  midpoint  of  the  breadth  B  of  the  cylinder  at 
the  calm-water  line.  The  cylinder  contour  below  the  calm-water  line  in  the  ly-plane  is  denoted 
by  S.  The  direction  of  propagation  of  the  incident  wave,  the  sinusoidal  wave  of  frequency  w, 
is  at  an  angle  ^  to  the  z-axis. 

In  the  linear  diffraction  problem,  which  WADAl  solves,  the  total  wave  potential  $vv  is 
separated  into  the  incident-wave  potential  $/  and  the  diflfracted-wave  potential 


^w{x,y,t)  =  +  ^D{x,y,t). 

(1) 

The  incident-wave  potential  is  given  by 

^/(a:,y,0  =  <pi{x,y)e~*'^^ 

(2) 

with 

<pj{x,y)  = 

t. » 

(3) 

where  i  =  t  is  the  time,  a  is  the  amplitude  of  the  incident  wave,  g  is  the  gravitational 

acceleration,  and  u  =  w*/y  is  the  wave  number.  Here  and  in  the  following,  it  is  to  be  understood 
that  only  the  real  parts  of  complex  expressions  are  to  be  taken. 
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The  diftracLed-wave  potential  is  assumed  to  be  of  the  form, 


(x,  y,  t)  =  <pd  (x,  y)e 


(4) 


where  fpD  is  the  solution  of  the  Laplace  equation: 

d^<PDix,y)  ,  d'‘<pDix,y)  ^ 

- ;r-5 - 1 - I - =0  in  the  fluid  domain,  y  <  0  (5) 

ax^  ay^ 

and  satisfies  the  following  boundary  conditions: 


d<PD 

dy 


VifD  =  0 


lim  { 

*— *±oo 


lim  — —  =  0 
v-^-oo  oy 

d‘PD  ^  - 

T  ii^<PD}  =  0 


on  the  free  surface,  y  =  0 

on  the  sea  bottom 

at  infinity  (radiation  condition) 


d(pD  _  d(p[ 
dn  dn 


on  the  body  surface  S(i,y)  —  0. 


(6) 

(7) 

(8) 
(9) 


In  Equation  (9),  djdu  =  n  •  V  denotes  the  derivative  in  the  direction  of  n,  the  vector 
normal  to  the  surface  and  pointing  into  the  body.  Except  for  the  body-boundary  condition 
(9),  the  above  boundary- value  problem  is  identical  to  the  strip-theoretical  radiation  problem 
in  which  the  right-hand  side  of  Equation  (9)  is  replaced  by  the  prescribed  velocity  of  the  body 
surface. 

In  oblique  waves,  <pd  should  depend  on  the  z-coordinate  as  well  as  on  the  x-  and  y- 
coordinates,  but  the  dependence  olipu  on  the  z-coordinate  is  neglected  in  the  two-dimensional 
mathematical  model  (5)-(9).  As  a  result,  the  diffraction  pressure  predicted  by  WADAl  is 
expected  to  become  less  reliable  as  the  direction  of  propagation  of  the  incident  wave  deviates 
away  from  that  of  the  beam  wave  towards  those  of  the  head  or  following  waves,  that  is, 
1/3  -  Jr/2|  ->  ir/2. 


3.  SOLUTION  TECHNIQUE 

A  common  scheme  of  solving  the  radiation  and  diffraction  problems  is  the  source-distribu¬ 
tion  method  using  the  Green  function  G{x,y\^,rj)  that  defines  a  velocity  potential  at  a  field 
point  (x,y)  owing  to  a  source  of  unit  strength  located  at  a  point  If  the  strength  of 

the  source  is  Q(^,r;),  then  by  the  assumption  of  linearity,  the  velocity  potential  at  (i,y)  is 
G{x,y;^,r))Q{^,r)).  Using  a  Green  function,  the  diffracted-wave  potential  <pD{x,y)  can  be 
expressed  in  terms  of  a  distribution  of  sources  of  strength  on  the  immersed  surface  5 

of  the  body  under  the  calm  water  line;  that  is, 

<PD{x,y)  =  G{x,y,^,r))Q{^,r))ds{^,r)).  (10) 

Implicit  in  Equation  (10)  is  the  dependence  of  <pd  on  i/  and  /3. 
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The  Green  function  G  used  in  WADAl  is  a  particular  solution  of  the  Laplace  equation, 
which  is  regular  throughout  the  fluid  domain  except  at  the  source  point  and  directly 

satisfies  the  boundary  conditions  on  the  free  surface  (6),  on  the  sea  bottom  (7),  and  at  infinity 
(8).  Denoting  a  field  point  by  z  ~  x  +  iy  (not  to  be  confused  with  the  z-coordinate  in  Fig.  1) 
and  a  source  point  by  f  =  ^  +  ir?,  G  can  be  written  as  (see  Reference  2); 


iog(z  -  f )  -  -  ?)  +  2P.V.  y 


roo 


1/  -  k 


■dk 


—  »Re 


(11) 


where  Re  indicates  the  real  part,  f  is  the  complex  conjuate  of  f  and  P.V.  indicates  the  Cauchy 
principal  value. 

The  unknown  source  distribution  along  S  is  determined  by  applying  the  kinematic 

condition  on  the  body  surface  (9).  Substituting  Equation  (10)  into  Equation  (9)  yields  the 
following  integral  equation  for  Q: 


ig(i,t/)4-P.v.|  QU,n) 


dG(x,y;^,rj)^^^^  ^  _^d^/(x,y) 


dn{x,y) 


dn(x,y) 


(12) 


Thus,  the  use  of  the  Green  function  (11)  reduces  the  boundary- value  problem  (5)-(9)  in  an 
infinite  exterior  domain  to  one  involving  only  the  body  surface  S,  Equation  (12). 

For  an  arbitrary-shaped  body,  the  integral  equation  (12)  can  be  solved  only  through 
a  numerical  method.  WADAl  uses  a  method  of  discretization  known  as  Frank’s  close-fit 
method  (Reference  1)  in  which  a  section  contour  5  is  represented  by  JV  straight-line  segments 
Sj(j  =  1,2,...,A^),  and  the  source  strength  Qj  along  an  individual  segment  is  assumed  to  be 
constant.  The  nodes,  or  the  points  at  which  the  unknown  values  are  sought,  are  taken  to  be 
at  the  midpoint  of  each  segment. 

Let  denote  the  complex  source  strength  of  the  segment  in  the  mode  of 

excitation,  where  m  =  1,2  for  the  horizontal  excitation  (asymmetric  abcut  the  y-axis)  and 
vertical  excitation  (symmetric  about  the  y-axis)  respectively.  Then,  for  every  and  /3,  the 
integral  equation  (12)  can  be  reduced  to  a  set  of  N  linear  algebraic  equations; 


(13) 


where  are  the  influence  coefficients  defined  by. 


(14) 


with  Sij  denoting  the  Kronecker  delta  function  {Sij  =  1  if  i  —  j  and  =  0  otherwise)  and 

Ktj  ^  u^-\/G{x„y,,^,r})ds{^,r]).  (15) 

The  components  of  the  normal  velocity  in  the  mode  of  the  incident  wave  at  5, 

are  found  from  (3): 


C//')  =  -n.  •  v(Re^/) 


(^.  ,v. ) 


=  [sin /9cos(i/a:,  sin /9)n,i  +  sin(p'a:,  sin  ,3)n,2] 


(16a) 
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=  -tn,  •  v(Im^;) 


(li.Vi) 


=  -towe*'*'’ [sin  ^ sin(i/i, sin  ^)n,i  -  cos(i/ri  sin  ^)n,2]- 


(166) 


Here  (li.yi)  are  the  coordinates  of  the  node  situated  at  the  midpoint  of  the  segment  Si, 
and  n»  =  (n,i,n,2)  is  the  normal  unit  vector  to  5,  pointing  out  of  the  fluid.  For  the  details  of 
evaluation  of  the  influence  coefficients  lij,  see  Reference  1. 

By  solving  the  algebraic  equation  system  (13),  the  source  density  can  be  deter¬ 

mined.  Then  the  diffracted-wave  potential  (pD{^t^y^)  for  is  found  by  a  direct  quadrature 
via  Equation  (10),  so  that 


where 
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FD(a:.,y.)  ^ 

Wl  1 

N  2 

=  DE‘?S">)r., 

y=l  m=l 

(17) 

r.j  =  G{x„yi,^,t})ds{^,r)). 

(18) 

The  total  periodic  wave-induced  pressure  Pw  is  obtained  from  the  linearized  Bernoulli 

equation,  Pw  =  -p~z —  for  each  node  (i«,y,),  with  given  in  Equation  (l).  Thus, 
at 


(it , yt , ()  —  Pi "t"  Poi^ityt!^) 

-  b/(a:.,y.)  +PD(^t,y,)je~“^‘ 


(19) 


where  p;  and  po — the  complex  amplitudes  of  the  incident-wave  and  diffracted-wave  pressures, 
respectively — are  given  by 


P/(a:.,y.)  =  iwpv^/(i,,y,) 
PD{x%,y,)-  «wp¥?r)(i,,y,) 


(20) 

(21) 


and  [Pwl  and  4>  are  respectively  the  magnitude  and  argument,  or  the  amplitude  and  phase 
angle,  of  Pw  ■ 


4.  IRREGULAR  FREQUENCIES 

Although  physical  problems  always  have  unique  solutions,  the  corresponding  integral  equa¬ 
tions  do  not  have  this  property  and  break  down  at  a  discrete  infinite  set  of  wave  frequencies 
called  irregular  frequencies.  This  phenomenon  was  first  pointed  out  by  John®.  To  determine 
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the  irregular  frequencies  for  an  arbitrary-shaped  body,  it  is  necessary  to  solve  the  adjoint  in¬ 
terior  problem;  see,  for  example.  References  1,  3,  and  4.  (For  a  rectangular  cylinder  such 
as  shown  in  Fig.  3,  the  irregular  frequencies  can  be  expressed  in  a  closed  form  as  follows^: 

=  sj  {nng/ B)  coth(nxT  /  B),  where  is  the  irregular  frequency,  n  =  1,2,3,....  Here, 
B  is  the  breadth  and  T  is  the  draft  of  the  cylinder.) 

In  practice,  the  influence-coefficient  matrix  (14),  which  approximatec  Equation  (12)  at 
discrete  points,  is  often  ill-conditioned  even  when  u  lies  in  the  vicinity  of  the  irregular  fre¬ 
quencies.  To  guard  against  the  occurrence  of  this  phenomenon,  the  user  should  carry  out  the 
computation  at  a  few  extra  frequencies  in  the  neighborhood  of  the  frequencies  of  interest  to 
make  sure  that  no  sudden  discontinuities  in  the  computed  values  occur  when  they  are  plotted 
against  frequency.  The  user  should  disregard  all  results  related  to  the  discontinuities  if  they 
should  occur. 


5.  DESCRIPTION  OF  INPUT 


Data  Set  1 
TITLE 

Data  Set  2 
ISYM 


NP 


Data  Set  3 

(XA(I),  1=  1,NP) 


Data  Set  4 

(YA(I),I=  1,NP) 


(72  characters  or  less  -  FREE  FORMAT) 

-  Brief  description  for  identification  purposes. 

(2  integers  -  FREE  FORMAT) 

-  Control  integer  indicating  whether  the  cross  section 
possesses  a  symmetry  with  respect  to  the  yr-plane: 

ISYM  =  0  No  left-right  symmetry; 

ISYM  =  1  Left-right  symmetry  exists. 

-  Number  of  offsets. 

Note; 

The  number  of  segments  should  never  be  less  than  six  on 
a  half-section  contour  to  ensure  sufficient  accuracy.  (The 
size  of  a  segment  must  be  much  less  than  the  wavelength 
and  the  local  radius  of  curvature  of  the  sectional  contour. 
Also,  because  the  source  strengths  are  averaged  over  each 
segment,  the  segment  sizes  should  be  made  smaller  near 
sharp  corners.) 

(NP  reals  -  FREE  FORMAT) 

-  i-coordinates  of  the  offsets. 

(Note  the  remarks  described  in  Data  Set  4  below) 

(NP  reals  -  FREE  FORMAT) 

-  y-coordinates  of  the  offsets. 
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Note: 


(1)  If  ISYM  =  0,  then  the  first  offset  (XA(l),YA(l))  is 
the  intersection  of  the  free  surface  and  the  body  contour 
in  the  left  half-section,  and  going  along  in  a  counter¬ 
clockwise  direction,  the  last  offset  (XA(NP),YA(NP))  is 
the  intersection  of  the  calm-water  line  and  the  body 
contour  in  the  right  half-section.  The  intersection  of  the 
contour  and  the  j/-axis  must  be  included  as  an  offset. 

(2)  If  ISYM  =  1,  then  the  first  offset  (XA(l),YA(l))  is  the 
intersection  of  the  y-axis  and  the  body  contour,  and  going 
along  in  a  counter-clockwise  direction,  the  laist  offset 
(XA(NP),YA(NP))  is  the  intersection  of  the  calm-water 
line  and  the  body  contour  in  the  right  half  section. 

(3)  The  computational  algorithm  in  WADAl  requires  that 
no  two  offsets  in  a  half-section  have  identical  vertical  co¬ 
ordinates.  For  horizontal  segments,  this  difficulty  can  be 
easily  overcome  by  giving  insignificantly  small  “artificial 
slopes”  to  the  segments;  compare  Fig.  3(a)  and  Appendix 
A(b). 


Data  Set  5 
NOX 


(1  integer  -  FREE  FORMAT) 
-  Number  of  wave  frequencies. 


Data  Set  6 


(NOX  reals  -  FREE  FORMAT) 


(XIB(I),  1=  l,NOX) 


-  Values  of  the  nondimensional  frequency  parameter: 

=  — — ,  where  w  is  wave  frequency,  g  the  gravitation- 
2^ 

al  acceleration,  and  B  the  waterline  breadth  of  the  body. 


Data  Set  7 


(1  integer  -  FREE  FORMAT) 


NOB  -  Number  of  wave  directions. 

Data  Set  8  (NOB  reals  -  FREE  FORMAT) 

(BETA(I),  1=  l,NOB)  -  Directions  P  of  wave  propagation  in  degrees;  see  Fig.  1. 


Sample  input  is  included  in  Appendix  A. 


6.  DESCRIPTION  OF  OUTPUT 

Appendix  B  shows  an  example  of  output  corresponding  to  the  sample  input  in  Appendix 
A.  The  first  part  of  the  output  is  a  printout  of  the  input  variables.  They  are  followed  by  the 
values  calculated  by  WADAl  for  every  combination  of  the  wave  frequency  parameter  (XIB) 
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and  the  wave  direction  0  (BETA).  The  real  and  imaginary  parts  of  the  complex  amplitudes 
of  the  incident-wave  pressure  pj  (PI)  and  the  diffracted-wave  pressure  pr>  (PD)  are  printed, 
followed  by  the  magnitude  (MAG)  \Pw  \  and  the  argument  (ARG)  4>  of  the  total  wave-induced 
pressure  Pw  (PW)  for  each  contour  segment.  Figures  2  and  3  show  the  plots  of  |Piv|  and  4>. 

7.  CONCLUDING  REMARKS 

The  FORTRAN  program  WADAl  was  developed  for  use  in  fluid-structure  interaction  prob¬ 
lems.  Using  the  source-distribution  method,  WADAl  calculates  the  hydrodynamic  pressure 
exerted  by  regular  incident  waves  and  diffreicted  waves  on  arbitrary-shaped  two-dimensional 
bodies  fixed  in  the  free  suface  of  deep  water.  WADAl  requires  a  simple  data  input  and  in¬ 
volves  a  small  number  of  unknowns.  This  is  because  the  Green  function  used  in  WADAl  is 
a  particuleir  fundamental  solution  that  directly  satisfies  the  boundary  conditions  on  the  free 
surface,  sea  bottom,  and  at  infinity,  leaving  only  the  body  boundary  to  be  discretized. 

Since  the  solution’s  dependence  on  the  longitudinal  coordinate  is  neglected  in  the  two- 
dimensional  mathematical  model  for  WADAl,  the  reliability  of  WADAl ’s  prediction  is  ex¬ 
pected  to  deteriorate  as  the  direction  of  propagation  of  the  incident  wave  deviates  away  from 
the  direction  of  the  i-axis. 
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FIG.  1.  Body  geometry  and  definition  of  coordinates 
(a)  waterline  plane  and  (b)  cross  section. 
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FIG.  2.  Results  for  an  asymmtric  cylinder  for  =  10>  ^  =  90°:  (a)  polygonal 
approximation  of  contour  and  offsets,  (b)  calculated  pressure  amplitude  around 
the  contour,  and  (c)  pressure  phase. 


0 (DEGREES) 

(b) 


FIG.  3.  Results  for  a  symmtric  cylinder  for  =  1-0)  —  135°:  (a)  polygonal 

approximation  of  contour  and  offsets,  (b)  calculated  pressure  amplitude  around 
the  contour,  and  (c)  pressure  pheise. 
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APPENDIX  A:  SAMPLE  INPUT 
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Data 

Set 

1 

FLOATING  ASn**(rrRIC  CYLINDER 

Data 

Set 

2 

0  17 

Data 

Set 

3 

-1.  -0.9  -0.8  -0.7  -0.5  -0.3  -0.15  0. 

0.1  0.2  0.4  0.55  0.6  0.7  0.8  0.9  1.0 

Data 

Set 

4 

0  -0.1  -0.2  -0.4  -0.58  -0.7  -0.72  -0.78 
-0.79  -0.8  -0.72  -0.60  -0.45  -0.30  -0.20  -0 

Data 

Set 

5 

1 

Data 

Set 

6 

1.0 

Data 

Set 

7 

1 

Data 

Set 

8 

90. 

(a)  AN  ASYMMETRIC  CYLINDER  (SEE  FIG.  2) 


Data 

Set 

1: 

FLOATING  RECTANGULAR  CYLINDER 

Data 

Set 

2: 

1  16 

Data 

Set 

3: 

0.  .156  .313  .469  .625  .781  .938  1.094  1.25 
1.25  1.25  1.25  1.25 

Data 

Set 

4: 

-1.001  -1.000875  -1.00075  -1.000625  -1.0005 
-1.000125-1.  -0.857  -0.714  -0.571  -0.429 

Data 

Set 

5; 

1 

Data 

Set 

6: 

1 . 

Data 

Set 

7: 

1 

Data 

Set 

8: 

135. 

-«.286 


1  .25 

-1 .00025 
-0.143  0. 


(b)  A  SYMMETRIC  CYLINDER  (SEE  FIG.  3) 


APPENDIX  B:  SAMPLE  OUTPUT 
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calculation  of  wave-induced  pressures  on  a  cylinder  fixed  in  the  free  surface 


INPUT  VALUES 


FILEIN  -  WADA1.CYL  IP  -  1 
TITL  -  FLOATING  ASYVWETRIC  CYLINDER 


SYM  -  0 

NOP  -17 

OFFSET 

XA 

YA 

1 

-1 .00000 

0 . 00000 

2 

-0.90000 

-0.10000 

3 

-0 . 80000 

-0.20000 

4 

-0.70000 

-0 . 40000 

5 

-0 . 50000 

-0 . 58000 

6 

-0.30000 

-0.70000 

7 

-0.15000 

-0 . 72000 

8 

0 . 00000 

-0.78000 

9 

0.10000 

-0.79000 

10 

0.20000 

-0 . 80000 

11 

0 . 40000 

-0.72000 

12 

0.55000 

-0 . 60000 

n 

0 . 60000 

-0 . 45000 

14 

0 . 70000 

-0.30000 

15 

0 . 80000 

-0.20000 

16 

0.90000 

-0 . 1 0000 

17 

1 . 00000 

0 . 00000 

NOX  -  1 

XI8  VALUE(s)  -  i.cecee 
NOB  *  1 

BETA  VALUE(S)  -  90.0 


CALCUUTED  VALUES 


HALF  BREADTH  -  1.00000 


PRESSURE  DISTRIBUTION  XIB  -  1.00000  BETA  -  90.0 


SEG 

REAL(PI) 

IMAG(PI) 

REAL(PD) 

IMAG(PD) 

MAG(PW) 

ARG(PW) 

1 

0.5533 

-0.7737 

-0.3285 

-0.7097 

0.2338 

15.8948 

2 

0.5681 

-0.6466 

-0.3655 

-0.6131 

0.2053 

9.3931 

3 

0.5420 

-0.5050 

-0.3611 

-0.5032 

0.1809 

0.5587 

4 

0 . 5056 

-0.3459 

-0.3532 

-0.4177 

0.1685 

-25.2313 

5 

0.4857 

-0.2053 

-0.3484 

-0.3706 

0.2149 

-50.2804 

6 

0.4793 

-0.1097 

-0.3453 

-0.3669 

0.2900 

-62.4862 

7 

0.4710 

-0.0354 

-0.3353 

-0.3609 

0.3527 

-67.3640 

8 

0.4555 

0.0228 

-0.3282 

-0.3590 

0.4025 

-71 .5523 

9 

0.4465 

0.0675 

-0.3063 

-0.3986 

0.4867 

-73.2628 

10 

0.4468 

0.1382 

-0.2940 

-0.4235 

0.5621 

-74.7834 

11 

0.4596 

0.2364 

-0.2850 

-0.4851 

0.7423 

-76.3934 

12 

0.4964 

0.3217 

-0.2925 

-0.5765 

0.9211 

-77.2115 

13 

0.5471 

0.4159 

-0.3205 

-0.6086 

1.0493 

-77.5263 

14 

0.5698 

0.5309 

-0.3201 

-0.6211 

1.1787 

-77.7669 

15 

0.5681 

0.6466 

-0.2916 

-0.6487 

1.3245 

-77.9519 

16 

0.5533 

0.7737 

-0.2416 

-0.7037 

1.5100 

-78.0873 

(a)  ASYMMETRIC  CYLINDER  (SEE  FIG.  2) 
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CALCULATION  OF  WAVE-INDUCED  PRESSURES  ON  A  CYLINDER  FIXED  IN  THE  FREE  SURFACE 


INPUT  VALUES 


FI  LEIN  -  WADA1.RAT  IP  -  1 
TITL  -  FLOATING  RECTANGULAR  CYLINDER 


SYM  -  1  NOP  -16 


OFFSET 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


XA 

0.00000 
0.15600 
0.31300 
0.46900 
0.62500 
0.78100 
0.93800 
1 . 09400 

1 .25000 

1.25000 

1 .25000 

1 .25000 

1 .25000 

1 .25000 

1 . 25000 

1 . 25000 


YA 

-1.00100 
-1 .00087 
-1 .00075 
-1 .00062 
-1 .00050 
-1 .00038 
-1  .00025 
-1.00013 
-1 .00000 
-0.85700 
-0.71400 
-0.57100 
-0.42900 
-0.28600 
-0.14300 
0.00000 


NOX  -  1 

XIB  VALUE(S)  -  1.00000 


NOB  -  1 

BETA  VALUE(S)  -  135.0 


CALCULATED  VALUES 


HALF  BREADTH  -  1.25000 


PRESSURE  DISTRIBUTION  XIB-  1.00000  BETA  -  135.0 


SEG 

REAL(PI) 

IMAG(PI) 

REAL (PD) 

IMAG(PD) 

MAG(PW) 

ARG(PW) 

1 

0 . 4486 

0.0196 

-0.3853 

-0.1336 

0.1659 

-67.5952 

2 

0.4451 

0.0594 

-0.3673 

-0.1115 

0.1877 

-65.5182 

3 

0.4381 

0.0985 

-0.3415 

-0.0922 

0.2138 

-63.1248 

4 

0.4278 

0.1368 

-0.3074 

-0.0771 

0.2454 

-60 . 6304 

5 

0.4141 

0.1740 

-0.2640 

-0.0679 

0.2846 

-58.1790 

6 

0.3972 

0.2099 

-0.2093 

-0.0672 

0.3348 

-55.8626 

7 

0.3771 

0.2442 

-0.1393 

-0.0798 

0.4019 

-53.7336 

8 

0.3541 

0.2765 

-0.0406 

-0.1207 

0.5061 

-51.7184 

9 

0.3617 

0.3091 

0.1194 

-0.2582 

0.7438 

-49.7001 

10 

0.4055 

0.3465 

0.1925 

-0.3441 

0.9136 

-49.1112 

11 

0.4547 

0.3885 

0.2508 

-0.4176 

1.0712 

-48.8085 

12 

0.5096 

0.4355 

0.3050 

-0.4894 

1 .2324 

-48.6286 

13 

0.5711 

0.4880 

0.3581 

-0.5631 

1 . 4030 

-48.5206 

14 

0 . 6404 

0.5472 

0.4118 

-0.6405 

1 .5867 

-48.4607 

15 

0.7180 

0.6135 

0.4665 

-0.7223 

1.7854 

-48.4367 

(b)  SYMMETRIC  CYLINDER  (SEE  FIG.  3)' 


-1 

-2 

-3 


-5 

-6 

-7 

-a 

-9 

-10 

-11 

-12 

-13 

-14 

-15 


0.4486 

-0.0198 

0.4451 

-0.0594 

0.4381 

-0.0985 

0.4278 

-0.1368 

0.4141 

-0.1740 

0.3972 

-0.2099 

0.3771 

-0.2442 

0.3541 

-0.2765 

0.3617 

-0.3091 

0.4055 

-0.3465 

0.4547 

-0.3885 

0.5096 

-0.4355 

0.5711 

-0 . 4880 

0.6404 

-0.5472 

0.7180 

-0.6135 

-0.3960 

-0.1573 

-0.3996 

-0.1818 

-0.3963 

-0.2062 

-0.3861 

-0.2293 

-0.3689 

-0.2503 

-0.3443 

-0.2680 

-0.3114 

-0.2805 

-0.2663 

-0.2826 

-0.2232 

-0.2597 

-0.2320 

-0.2654 

-0.2493 

-0.2810 

-0.2719 

-0.3032 

-0.2997 

-0.3318 

-0.3328 

-0.3669 

-0.3715 

-0 . 4092 

0.1472  -69.0807 
0.1306  -69.6301 
0.1155  -68.7536 
0.1015  -65.7268 
0.0888  -59.3571 
0.0786  -47.7169 
0.0750  -28.8999 
0 . 0880  -3 . 9509 
0-1471  19.6045 
0.1916  25.0661 
0.2319  27.6386 
0.2720  29.0904 
0.3132  29.9311 
0.3566  30.3797 
0.4022  30.5362 


(b)  SYMMETRIC  CYLINDER  -  CONTINUED 
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